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The generalized Chaplygin gas, characterized by the equation of state p = —A/p a , has been 
considered as a model for dark energy due its dark-energy-like evolution at late times. When 
dissipative processes are taken account, within the framework of the standard Eckart theory of 
relativistic irreversible thermodynamics, cosmological analytical solutions are found. Using the 
truncated causal version of the Israel-Stewart formalism, a suitable model was constructed which 
• cross the w = — 1 barrier. The future-singularities encounter in both approaches are of a new type, 

not included in the classification presented by S. Nojiri, S. D. Odintsov and S. Tsujikawa, Phys. 
O ! Rev. D71, 063004 (2005). 
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I. INTRODUCTION 

£C) ' In the framework of general relativity, the acceleration in the expansion of the universe during recent cosmological 
times, first indicated by Supernovae observations [l| and also supported by the astrophysical data obtained from 
WMAP indicates the existence of a exotic fluid with negative pressure, which constitutes about the 70 per cent of the 
total energy of the universe. To know the nature and the behavior of this dark energy is one of the great challenges 
of cosmology and, indeed, of fundamental physics. 

Within the different candidates to play the role of the dark energy, the generalized Chaplygin gas (GCG), in which 
^ i the pressure p and the energy density p are related by 

p = -A/ P a , (1) 

bJ)' where A and a are constants, has emerged as a possible unification of dark matter and dark energy, since its cosmo- 
logical evolution is similar to an initial dust like matter and a cosmological constant for late times. This idea was first 
proposed inQ • Using supernova observations this model has showed to be consistent for any value of a in the range 
< a < 1 0f 

Nevertheless, L. Amendola et al jj| found that the latest CMB anisotropy data seems favoured GCG as a dark 
energy model instead of a candidate for a unified dark model. The parameter a is rather severely constrained, i.e., 
< a < 0.2 at the 95% confidence level. Others works, such as @ tends to confirm that adiabatic Chaplygin gas 
is ruled out, as unified model, unless the parameter a be very close to zero. A Chaplygin quartessence model with 
vanishing pressure perturbation, the called Silent Chaplygin, has shown to be consistent with CMB data Q . 

Since we are interested here only in the late time acceleration and the phantom behavior, we will considered a 
GCG as a model for dark energy. Specifically, our principal aim is to investigate the possibilities that GCG models 
could realize the phantom divide u> = — 1 crossing phenomenon. Previous works have indicated this possibility. For 
example, using an extension of the Chaplygin gas model, where A was taken as a function of the scale factor, a, of 
the form A/a~ m ( A/ is a constant and m > 0), Meng et al Q showed that exist some cases which can realize state 
parameter u = — 1 crossing. 

Other approaches have considered dissipative effects in GCG models, using the framework of the non causal Eckart 
theory |e| ■ In Q a viscous GCG was investigated from the point of view of the dynamical analysis, assuming that 
there is bulk viscosity in the linear barotropic fluid and GCG. It is found that the equation of state of GCG can cross 
the boundary w = — 1. 
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It is interesting to note that the a GCG itself may behave like a fluid with viscosity, in the context of the Eckart 
formalism. J. C. Fabris et al [Tof found an equivalence between the GCG and a dust like fluid with an special election 
of both parameter m and a, where m is related with the viscosity coefficient throughout the relation = (3p m . 

It is well known that the Eckart 'approach is a non causal theory. The full causal theory developed by Israel and 
Stewart leads to stable behavior under a wide range of conditions and transient phenomena on the scale of the 
mean free path/time |l2j |. In this theory the character of the evolution equation is very complicated, nevertheless we 
were able to find solutions in the truncated case for dissipative GCG cosmological models. 

In the present paper we have investigated the behavior of a flat universe filled with a viscous GCG, using first the 
non causal approach of Eckart and then the truncated version of the Israel-Stewart formalism. In the context of the 
Eckart formalism, we show analytical solutions, which presents future-singularities, which were found restricting the 
range of the parameters m and a. 

Using a truncated version of the Israel-Stewart formalism we were able to find a suitable model that present 
future-singularities, but with an equation of state corresponding to quintessence. This model has a relaxation time 
compatible with the expansion of the universe in a regime near to the thermodynamic equilibrium. 

The organization of the paper is as follows: in Section II we present the field equations for a flat FRW universe 
filled with a viscous GCG within the framework of the Eckart theory. Cosmological solutions are obtained solving a 
non linear differential equation for the Hubble parameter. Exact solutions are obtained in terms of LerchPhi function 
for m and a satisfying the constraint given by m = — (a + 1/2). In Section III we solve the equations of the truncated 
version of the Israel-Stewart formalism, using an Ansatz for the Hubble parameter. In section IV we characterize the 
future-singularities obtained. In section V we summarize our results. 



II. THE GENERALIZED CHAPLYGIN GAS AND THE ECKART THEORY 



The FRW metric for an homogeneous and isotropic flat universe is given by 

ds 2 = -dt 2 + a(t) 2 (dx 2 + dy 2 + dz 2 ) , (2) 

where a(t) is the scale factor and t represents the cosmic time. In the following we use the units 8irG = 1. The field 
equations in the presence of bulk viscous stresses are 

-Y = = (3) 



- = H + H 2 = -Up + 3(P + Tl)), (4) 
a 6 

where n is the bulk dissipative pressure. The conservation equation is given by 

p + 3H(p + p + n) = 0. (5) 
Assuming that the dark component is GCG, the equation of state can be rewritten in the form 

p(p)=-M 4(1+c V Q , (6) 

where M has dimension of mass and a is a constant. The original Chaplygin gas correspond to the election a = 1. 
We note that a = — 1 yields the state equation corresponding to a cosmological constant. 

If we put the Eq. 10 in the barotropic form p(p) = up, the parameter u take the following form in terms of the 
Hubble parameter 

M 4 
3H< 



= " I ^) (7) 



From Eqs. © to JJJ) is direct to obtain a single evolution equation for H: 

2H + 3{l + Lu{H,a))H 2 = -U. (8) 
In the first order thermodynamic theory of Eckart 8] n is given by 

H = -3H£, (9) 
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and, in order to obtain solutions of the equation JSJ, we will assume that the viscosity has a power-law dependence 
upon the density 



£ = f3p m , /3>0. 



(10) 



where (3 and m are constant parameters. Using Eq. (JJJ), (|10J) and © in Eq. (JHJ), we obtain the non linear equation 
leading the evolution of the Hubble parameter 



2H + 3 <^ H 2 



M 4 \ 1+Q 

— J H- 2a + T l l3H 2m+1 



0. 



(11) 



Despite the complexity of this equation it is possible to obtain analytic solutions for special cases in which the 
parameters a, (3 and m are constrained. In the following subsections these solutions are explicitly showed. 



A. Solutions with m = — (a+ 1/2) 

If in Eq. (|ll|l we choose the exponents of H to be equal, i.e., —2a — 2m + 1 we obtain 

2H + 3H 2 - (3(a)H- 2a = 0, 

where (3(a) is defined by 

/3(a) = 3 1/2 - a 

Integrating Eq. (|12|) we obtain an implicit relation for H as a function of the cosmic time t 

t 



v3 



JLl^LerchPhi ( Mff-^+a) x __}_ 
3(1 + a) V 3 2(1 + a) 



where the LerchPhi function is defined as follows 



LerchPhi(z, a, v) = 



-! (n + v) a ' 



:i2) 



(13) 



(14) 



(15) 



71=0 



It was mentioned in that a GCG (without viscosity) is equivalent to a dust with viscosity, if the parameters a 
and m satisfy the constraint m = — (a + 1/2). In our case this constraint was imposed in order to obtain analytic 
solutions of a GCG with viscosity. 



B. Solutions with future-singularities in the non causal approach 

It is easy to show graphically that above solutions, for positive values of a, represents universes with future- 
singularities. Before to show an analytical solution of Eq. (|11|) . let us discuss some interesting results related to a 
pure GCG without dissipation. From Eqs. and (jSJ), we obtain the evolution of the energy density p in terms of the 
scale factor a, or equivalently in terms of the redshift z = a^/a — 1. The expression for p(z) is given by (for a ^= — 1) 



Af 4 l + g(l + z) 



3(l+a) 



l/(l+a) 



(16) 



where g is a parameter that can be fitted by observational data. 

Using the above equation in © we can obtain an expression for u> in terms of the redshift 



w(z,a) = -(l+. 9 (l + z) 3(1+a) ) 



(17) 



From Eq. (|17f) is direct to see that for negative z (future cosmic time), we can obtain u) = — 1. The phantom case 
can be obtained if g < 0. A. Sen and R. J. Scherrer 01 explore the possible scenarios when the usual constraint on 
the parameters a and g are relaxed. They consider the possibility of take g < and a < — 1. Two cases of phantom 
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dark energy are obtained. One of them is a late phantom GCG model which gives lu < — 1 with uj decreasing with 
time. A future singularity occurs at finite value of the cosmic time. 

In our scheme it is possible to obtain phantom evolution in the special case, which is straightforward to integrate, 
when m = 1/2. The bulk viscosity coefficient £ takes the form 



£ = (3p 1/2 = V3/3H. 

Introducing Eq. ljT%|) in Eq. (JSJ) , we obtain 

2H + 3 (l - y/Ep) H 2 = -3w (H) H 2 
which with the election (3=1/ the above equation has the solution 

r4 \ 1 + a 



H(t)=H 



l + -(2a+l) 



M 4 
3H 2 



-i 1/(2q+1) 



Assuming that 2a + 1 



-5 < 0, the solution 120(1 can be rewritten as 
H(t) = H ^ t^) 1 " 



obtaining, from Eq. J2J, the explicit expression for the energy density , 



which blows up at the time t s , given by 



p(t) = 3H 2 



2 

ts ~ 3d 



it - 1 ) 



-2/5 



3_Hl 

M 4 



(l-5)/2 



(18) 
(19) 

(20) 

(21) 
(22) 
(23) 



The solution given in Eq. (1221 implies that + H 2 > 0, i.e., the acceleration is positive. An expression for a(t) can 
be easily obtained integrating Eq. (|2*2*|l 



a (t) — a,Q exp 



-Hots 



1 - 



(* ~ h) 



-(l+S)/5 



- 1 



(24) 



which at t — t a — > t s blows up, representing a solution with a future-singularity. In brief, a GCG with bulk viscosity 
contains big rip solutions if a < —1/2 and if the bulk viscosity coefficient takes the specific form £ = -Lp 1 / 2 . 

It is interesting to compare the above results with those obtained in |13| , where the model with a future singularity 
occurs when g < and a < — 1. Nevertheless, it is easy to show from Eq. (|17|) that if g < it is possible to find, 
in the past or in the future evolution of the universe (depending on the value of g) that the fluid behaves as a stiff 
matter. This result is not physically plausible as we know from the cosmological observations that ruled out stiff 
matter, which would be of the importance only in the very early universe (see, for example |l4jp. In our case, a CGC 
fluid with bulk viscosity, the solution with future singularity is obtained for g > and a < —1/2. 

Note the exponential behavior of the scale factor with time. In the case where dark energy component obey the 
state equation p = (7 — \)p and £ ~ p 1 / 2 there are also big rip solutions but with a power-law expansion for the scale 
factor of the form 



a(t) = a ( 1 - 



t — to \ 3(T-V^i3) 

^7 



(25) 



where do = a(t = t 
and 7 are imposed |1 



This solution was found by Barrow |15| and if the following constrains on the parameters (3 

(26) 



V3/3> 7 , 

the scale factor blows up to infinity at a finite time t s > to, given by 

2 



ts 



3(V3/3 - 7) 



Ho- 1 - 



(27) 
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III. A TRUNCATED VERSION OF THE ISRAEL-STEWART FORMALISM 

In this approach we consider an imperfect fluid whose dissipative bulk pressure, II, obeys the causal transport 
Eq. [13 

n + rll = -3iJ£, (28) 

where r is the bulk relaxation time and is the crucial thermodynamic parameter in the Israel-Stewart formalism that 
ensures causality. 

In order to obtain an equation for II in terms of Hubble parameter H, we derive Eq. (JHJ) introducing then Eq. J2HJ- 
This yields 

11(H) = r (2H + 3uj(H, a)H 2 + 6(1 + u(H, a))HH^j - 3£H. (29) 
The parameter r is related to the energy density throughout the equation 

r=^=f3p m -\ (30) 
P 

where £ is given by Eq. {TDJ. In 

big rip solutions were found within the framework of the full causal theory of 
Israel- Stewart, for universes filled a dark component obeying the state equation p = (7 — l)p. For this case, using the 
Ansatz H = A(t s — f) , where A = Hot s and t s is the time when the future-singularity occurs, we were able to find 
the exact expression for A. 



A. Solutions with future-singularities in the causal approach 

In the following we will generalize the expression for the Ansatz thinking in the more complex structure of the 
equation of state corresponding to a GCG. Our election is 

H = AA~ q , (31) 

where A = H (t s — t ) q , A = t s — t and q is a positive constant, which is determined by consistence of the equations. 
Using the above Ansatz and Eqs. ©, 10 and (|TU|l . in Eq. we obtain 

n(A) = ^ f 3A- 2 ^ m+1 / 2 '> (a(m) - b(m, q)A q ~ 1 - c(m, q)A 2 ^ - ad(m, q; a)A^ 3+2a ^ , (32) 

where the coefficients a(m) and b(m, q) are given by the following expressions 

a{m) = 3 A 2 ( m+1 /2) ; 5( mj(? ) = 2 qA 2m , (33) 

and c(m, q) and d(m, q; a) by 

c(m,q) = -q{q + l)A 2 ^- l / 2 \ d(m, q; a) = 2q( ^- J ^2(m-x-a) < (34) 

Condition IT < 0. Since from thermodynamics arguments (as requiered from the second law), II < 0, we will explore 
the possible scenarios which are obtained depending on the values of the parameters q, m, and a. 

Case q = 1. Let us investigate the special case with q = 1, where the factors with the time dependence A 9-1 in 
Eq. (|32|l . becomes constant. In this situation we obtain constraint for A which are independent of the parameter m. 
In brief, we have the following possibilities in order to ensure II < 0: 

A > a = 0, (35) 



and 



A> a<0 - ( 36 ) 



6 



In the first case it is straightforward to see that the time t s has a lower limit given by 

t. > ^W 1 - (37) 

Case q > 1. Choosing for simplicity a = 0, the sign of Eq. Ij32|l becomes dependent of the sign of the factor 
a(m) — b(m,q)A q ~ 1 — c(m, q)A 2 ( g ~ 1 ' 1 . Since the powers of the A functions are both positives, they are decreasing 
functions of the cosmic time. This implies that if we demand the above factor to be positive at t = it will positive 
for the future evolution. Solving the constraint 

(a(m) - b(m, q)A q - 1 - c(m, q)A 2 ^) | t=0 > 0, (38) 

and since a(m),b(m,q),c(m,q) are functions of A — Ho(t a — to) q , at t = Eq. I|38(l gives a constraint for t a in the 
form of V(t s ) > 0, where V(t s ) is a polinomial of the form 

Pxtl-fatlT 1 -fot-\ (39) 

which can be solved for t s . 

Case < q < 1. For this case II < if a < 

Condition |II| << \p\. We have just found the constraints on the parameters a, depending of which values are 
given for q, in order to obtain a negative dissipative pressure. Note that these constraints are independent of m. 
Constraints on m appears when we consider that the Israel-Stewart theory is derived under the assumption that the 
thermodynamical state of the fluid is close to equilibrium, which means that the nonequilibrium bulk viscous pressure 
should be small when compared to the local equilibrium pressure that is |I1| << \p\. Using Eqs. © and Q) the 
above condition takes the form 

l n ( A )l _ sm-lg f Ml\ ( + } J^a^-2q( m +l/2+a) x 



b(A)| 1 V 3 

(a(m) - b{m, q)A q - 1 - c(m, q)A 2 ^-^ - ad(m, q; a)A q( - 3+2a ^ 1 ^ « 1. (40) 

It is straightforward to check out that both conditions, II < and |IT| << \p\ are satisfied in the following cases: 
Case q = 1 



and 



a = 0, m+i<0, (41) 



A = l - ± ^ L i «<« ™<\- ( 42 ) 



Case q > 1 



a = 0, m+-<0. (43) 

We finally found the value of the parameter (3, constraining the relaxation time r to be lower than the Hubble time, 
i.e., t < H . The expression for r is given by 

t = ~ = 3 m - 1 / g J 4 2 ( m - 1 / 2 )/\- 2 9( m - 1 / 2 )/f- 1 (44) 
P 

Since the relation r < H -1 must be satisfied for all times, at t = we obtain and upper limit for the parameter /3, 
which can be evaluated from 

(S'"- 1 ^ 2 ^- 1 / 2 ^- 2 ^" 1 - 1 / 2 )) | t=0 < 1, (45) 

in terms of t s . In order to satisfy r < H^ 1 for all times we need m < 1/2. If this condition holds, r is a decreasing 
function of the cosmic time and becomes equal to zero at t s . We can take, for example, a = 0, m < — 1/2 (from the 
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condition (|41|l . In other words, with m < —1/2 (and q = 1 and a = 0), it is possible to have a realistic models with 
IT < 0, |II| << p and r < H^ 1 for all times during the cosmic evolution. 

Let us show the behaviour of one of these models, taking for simplicity, q = 1, a = and m = — 1. In this model 
the solution for the Hubble parameter is given by 

H = AA- 1 , (46) 
From Eq. we obtain the expression for the dissipative pressure 

n(A) = - ^ A ~ A ^- A -\ , (47) 

where A± = 1/3(1 ± VE). As we saw above, A must satisfy the inequality given in (|35[l in order to ensure II < 0. 
The expression for pressure of the GCG in terms of used Ansatz is 

p(A) = -3f— J A -2<x A 2qa^ ( 4g ) 

which gives for the particular values that we are taking 

p(A) = -M A . (49) 

The negative pressure of the GCG will be constant during the entire cosmic evolution. The energy density can be 
obtained from Eqs. Q and lj4^|l . yielding 

p(A) = 3A 2 A" 2 , (50) 

which becomes infinity at t s . The upper limit that the parameter /3 could take can be evaluated from the conditions 
|II|/M << 1 and t < if -1 . Evaluating |II|/|p| << 1 for our model we obtain that 

|H(A)| _ (3 (A~A + )(A-A_) 

which implies the following upper limit for f3 when equation is evaluated at t = 

3/4 3 M 4 

^« { A-A +){ A-A A rni lH ° > - ' <52) 
The constraint r < H~ Y evaluated at t = gives the other restriction on f3 

f3 <9H$. (53) 

If we choose j3 constrained by Eq. H53JI . Eq. (|52J) can be easily satisfied. This means that with a enough low viscosity, 
the dark energy, modeled by a GCG, leads to a superaccelerated universe. 



IV. CHARACTERIZATION OF THE FUTURE-SINGULARITIES 



We will discuss briefly the conditions which presents the future-singularities founded in the above sections. In the 
two solutions found in section II, for a > 0, the future-singularities have the following characterization: 

• For t — > t s , a — > oo, p — > oo, and |p| — > 

The model explicitly calculated for the special values q = X,a = and m = —1, in section III, gives a solution 
which presents a future-singularity characterized by 

• For t — > t Sl a — ► oo, p — * oo, and \p\ — > constant 

In the classification realized in exist four types of singularities and none of them behaves like those founded 
in this work. In this sense, cosmological models filled with a dissipative GCG, in the framework of the truncated 
Israel-Stewart formalism, present new types of future-singularities. 
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V. CONCLUSION 

In the present paper we have found cosmological solutions for a GCG with bulk viscous stresses, in the context 
of both thermodynamic formalism of Eckart and Israel-Stewart. We have assumed that the viscosity has a power- 
law dependence upon the density, i.e., £ = (3p m , where (3 and m are constant parameters (/3 > 0). Following the 
Eckart approach and in order to obtain solutions, we have derived the non linear differential equation which leads 
the evolution of the Hubble parameter. Analytical solutions are found constraining the parameters a and m by the 
relation m — — (a + 1/2) and taking particular values of f3. Solutions with future-singularities can be obtained taking 
the specific value m = 1/2. In this case, the scale factor grows exponentially, differing from the power law expansion 
of the scale factor when the universe is filled with a fluid obeying a barotropic equation of state and viscosity. As it 
was discussed in 0], phantom scenarios with big rip are obtained if a < — 1 and g < 0. Nevertheless, those obtained 
from models with a GCG and bulk viscosity avoid the fact that the cosmic fluid could behave like stiff matter in the 
past or future evolution of the universe. 

Using the causal but reduced version of the Israel-Stewart formalism we have found a suitable physically solution 
from the point of view the thermodynamics requirements expressed in the conditions r < H^ 1 and |II|/|p| << 1. We 
have explicitly showed a case in which is possible to cross the barrier w = — 1. 

The future-singularities encounter in our models, in the framework of the non causal and causal formalism, are of 
a new type, not included in the classification presented in |Ts|. 
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